Astronomy & Astrophysics manuscript no. ms 


©ESO 2012 


October 3, 2012 





Uncovering the magnetic environment of our solar system 

(Research Note) 



L. Ben-Jaffel 1 ' 2 and R. Ratkiewicz 3 ' 4 



(N 

o 

-i— > 
O 

O 

(N 
P? 

6 

-t— > 

d 



1 UPMC Univ Paris 06, UMR7095, Institut d'Astrophysique de Paris, F-75014 Paris, France, e-mail: bjaffel@iap . fr 

2 CNRS, UMR7095, Institut d'Astrophysique de Paris, F-75014, Paris, France 

3 Institute of Aviation, Warsaw, Poland 

4 Space Research Center, Polish Academy of Sciences, Warsaw, Poland, e-mail: roma@cbk . waw . pi 



Received May 21, 2012; accepted August 11, 2012 



ABSTRACT 



Since its formation 4.6 billion years ago, our solar system has most likely crossed numerous magnetized interstellar clouds and bubbles 
of different sizes and contents on its path through the Milky Way. Having a reference model for how the heliosphere and interstellar 
winds interact is critical for understanding our current Galactic environment, and it requires untangling the roles of two major actors: 
the time-variable solar wind and the local interstellar magnetic field. Numerical simulations predict a distortion of the heliosphere 
caused by both solar wind anisotropy and interstellar magnetic field orientation. However, model comparison to deep space probes' 
measurements led to contradictory reports by Voyager 1 and Voyager 2 of both several crossings of the solar wind's termination shock 
and of the strength of the local interstellar field, with values ranging from 1.8 to 5.7 pG. Here, we show that Voyager 1 & 2 plasma, 
fields, and Lyman-a sky background measurements, as well as space observations of high-energy particles of heliospheric origin, may 
all be explained by a rather weak interstellar field 2.2 ± 0.1 jiG pointing from Galactic coordinates (/, b) ~ (28, 52) ± 3° . For the 2000 
epoch Ulysses-based helium parameters assumed thus far, the interstellar bow shock must exist. By contrast, using the 2010 epoch 
IBEX-based He parameters and a stronger magnetic field leads to a plasma configuration that is not consistent with the Voyagers TS 
crossings. For the newly proposed interstellar He parameters, more simulations are required before one may determine whether the 
interstellar bow shock truly does disappear under those assumptions. 

Key words, solar neighborhood — local interstellar matter — ISM: magnetic fields — Sun: heliosphere — solar wind — 
magnetohydrodynamics(MHD) 
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1. Introduction 

The supersonic, expanding solar wind (SW) slows through the 
termination shock (TS) to adjust to the outer conditions and 
then carves a cavity into the local interstellar medium (LISM), 
which is called the heliosphere . While the asymmetry of 
the heliosphere caused by the local interstellar magnetic field 
(ISMF) was fir st predicted by the Newtonian approximation 
(Fahr et al. 1988I) . it was not confirmed until 10 years later by 
full 3D magnetohydrod ynamic (MHD) theoretical models (e.g., 
(Ratk iewicz et al. 19 98)). Since then, modelers quickly realized 
that the observation/modelling of such asymmetry may con- 
strain the direction of the ISMF and act as an interstellar mag- 
netic compass. Before MHD models were employed, the first 
attempt to uncover the ISMF vector was made using observa- 
tions of an asymmetry in the spatial distribution of the Lyman- 
a emission far from the Sun (iBen-Jaffel et al. 20001) . Assuming 
the observed Lyman-c glow of heliospheric origin, the simple 
Newtonian approximation showed that the ISMF vector points 
~ 40° away from the LISM flow direction with a strength of 
~ 1.8 fiG. These results were later confirmed by MHD 3D simu- 
lations, providing an independent estimation of the ISMF direc- 
tion pointing from Gala ctic coordinates (37 + 14°, 49.5 + 8.5°) 
(iRatkiewicz et al. 20081) . 

Apart from the Lyman-a radiation constraints on the he- 
liosphere models, it is commonly believed that Voyager 1 
(VI) crossed the TS for the first time at t he heliocentric dis- 
tance ~ 94 AU on 2004 December 16 dBurlaga et al. 20051: 



iDecker et al. 20051: iGurnett & Kurth 20051: IStone et al. 20051) . 
and Voyager 2 (V2) at ~ 84 AU on 2007 August 30 
dBurlaga et al. 2008t IDecker et al. 200"1 IGurnett & Kurth 20081: 
iRichardson et al. 20081) . The 10 AU difference in the cross- 
ing distance indicates heliospheric asymmetries that could re- 
sult from the ISMF, the asymmetric SW dynamic pressure, or 
the motion of TS caused by any time-dependent phenomena 
(IStone et al. 20081) . 

More recently, IBEX has completed sky maps that im- 
age energetic neutral atoms (ENAs). These maps reveal a 
bright ribbon of ENAs in the energy range ~ 0.2 - 6 keV 
dMcComas et al. 20091 iFunsten et al. 20091 Fuselier et al. 20091: 
ISchwadron et al. 20091) . We note that Cassini observed 
~6- 1 3 keV ENAs with a similar but broader structure 
(iRrimigis et al. 20091) . One possible explanation assumes 
that the ribbon discovered by IBEX is pro bably ordered by the 
ISMF interacting with t h e heliosphere (McCo mas et al. 20091: 
iHeerikhuisen et al. 20101: IChalov et al. 201 OIK Because they 
provide global maps of the interstellar interaction, ENAs 
observations are highly complementary to and synergic with the 
detailed single in situ measurements provided by the Voyager 
probes. 

Most of these observations as well as further measurements 
of the plasma flow and magnetic field in the inner heliosheath 
by Voyager 2, have spurred numerous attempts to determine the 
parameters of ISMF from MHD and neutrals modeling. Table 
1 provides a list of papers in which the ISMF parameters have 
been estimated. A quick look shows that up to now no analysis 
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Table 1. Summary of models indicating the strength and orientation of the ISMF. 



ISMF strength QiG) 


ISMF orientation 


Be W 

U*, YY 


Solar-wind 


Lyman-O' 


VI 


V2 


IBEX 


Ref. 


1.8 


a ~ 40° 




SS 


+ 








(Ben-Jaffel et al. 2000) 




8 = 






-l- 








(Lallement et al. 2005) 


1.8 


(228 ± 14°, 42 ± 8.5°) 




ss 


+ 








(Ratkiewicz et al. 2008) 


3.8 


(248, 35)+5° 




SS 






+ 




(Ratkiewicz & Grygorczuk 2008) 


3.7-5.5 


a ~ 20 - 30, 5 < 90 


7 


ss 




+ 


+ 




(Qpher et al. 2009) 


> 4.0 


a ~ 30, /3 = 


+ 


ss 




+ 


+ 




(Pogorelov et al. 2009) 


3.0 


(224,41)° 


? 


ss 








+ 


(Heerikhuisen et al. 2010) 


4.4 


a ~ 20 


? 


ss 








+ 


(Chalov et al. 2010) 


3.0±1.0 


(225, 35)+5° 




ss 








+ 


(Grygorczuk et al. 2011) 


2.4+0.3 


(227, 35)+7° 


+ 


NS 






+ 


+ 


(Strumik et al. 2011) 


2-3 


(220-224, 39-44)° 


+ 


ss 


+ 






+ 


(Heerikhuisen & Pogorelov 2011) 


> 3.0 


a = 45° 


+ 


NS 








+ 


(McComas et al. 2012) 


2.2±0.1 


(224, 36)+3° 


+ 


NS 


+ 


+ 


+ 


+ 


This work 



Notes, a is angle between inflow V„ and magnetic field B„ vectors. /? is angle from HDP plane. For all references except McComas et al. (2012), 
the HDP was defined by Lallement et al. (2010). The ISM direction, speed, plasma temperature, and HDP orientation have different values in 
McComas et al. (2012). Coordinates are ecliptic. Bsw is the interplanetary magnetic field. 6 is the angle between the BV plane and the Sun's 
equatorial plane (noted ft in Opher et al. (2009)). Lyman-a corresponds either to Voyagers or to SOHO/SWAN ultraviolet observations. VI and 
V2 correspond to plasma and field in situ measurements. IBEX corresponds to the ENA ribbon observations. SS is for spherical symmetric solar 
wind flow and NS is for non-spherical symmetric. 



thus far proposed has used the full set of observations. In fact, 
with the last test reported, it was not possible to get the ISMF 
strength and orientation consistently using the crossing TS dis- 
tance from the Sun for VI ~ 94 AU, for V2 ~ 84 AU, and for the 
IBE X ribbon angula r location for the same initial S W conditions 
dStrumiket al. 201 lh . 

Here, we focus o n extending our previously reported work 
dStrumik et al. 20 lib , using a more accurate sampling of the pa- 
rameters space of our MHD model. Therefore, for the period of 
2002-2007, which covers both VI & V2 crossings, we report a 
rich set of numerical simulations that include several strengths 
and orientations of the ISMF, many SW conditions at differ- 
ent epochs that describe its anisotropy, and a fine 3D grid. Our 
primary goal is to fit VI, V2, and IBEX observations simulat- 
neously, which should provide an accurate measurement of the 
unperturbed ISMF vector at the actual Galactic position of our 
solar system. A comparison of our findings with the most recent 
reports on the ISMF vector is also provided. 

2. Model and method 

We use the 3D MHD time-dependent model in which the ef- 
fects of the interplanetary magnetic field and the anisotropic so- 
lar wind are simulated. We sol ve the set of MHD equations wit h 
a constant flux of hydrogen (Ratkiewicz and Ben-Jaffel 2002). 
This assumption is validated afterward by comparing our final 
results to independent model s that use a full kinetic descriptio n 
for neutrals (e.g., section 3) (Heerikhuisen & Pogor elov 201 lh . 
Furthermore, as a lever to uncover the ISMF vector, the TS 
asymmetry was shown to be unchanged by the treatment of neu- 
trals, provided the density level at that position is consistent wit h 
observations (lAlouani-Bibi et al. 201 \[ llzmodenov et al. 2 003). 
In our calculations the coordinate system is Sun-centered, and 
the x-y plane contains the ISMF B m and velocity V M vec- 
tors, with the latter parallel to the x-axis direction. At the 
outer boundary we use an LISM with plasma density n„ = 
0.095 cirT 3 and velocity V^, = 26.4 km/s, flowing from the 
He direction (255.4°, 5.2°) in ecliptic coordinates and tem- 
perature Too = 6400 K (her eafter called Ulysses-based He pa- 
rameters ( Witt e et al. 20041) ). To have the neutral H number 
density at the TS consistent with measurements, we choose 
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Fig. 1. Contours of regions in ecliptic coordinates of ISMF ori- 
entations where the MHD model reproduces the TS crossing for 
V2 on 2007 August 30 for ISMF strengths equal to (blue) 1.8, 
(green) 2.0, (orange) 2.2, (red) 2.4, and (yellow) 3 fiG (contours 
are shown from right to left in increasing strength). For the IBEX 
ribbon, the small black contour shows |B ■ R| < 0.025 obtained 
for 2.2yuG as the best fit. We note that the B ■ R = condition is 
not sensitive to the ISMF strength and only the best solution is 
shown. 



= 0.1 1 cm" 3 , which is still in the range 0.1 - 0.2 cm~ 3 indi- 
cated by (Izm odenov et al. 20031) for the LISM. The unperturbed 
ISMF strength varies within the limit 1.5 to 3.0 fiG. The orien- 
tation of the ISMF in space is defined by two angles: an inclina- 
tion angle a (between inflow V M and magnetic field B M vectors) 
which varies from 0° to 90°, and a deviation angle /? which varies 
from -90° to 90° and measures the an gle of the x-y plane f rom 
the hydrogen deflection plane (HDP) dLallement et al. 20101) . 

The anisotropy of the solar wind flow is included using slow 
and fast wind sectors that are estimated from yearly maps of 
observed interplanetary scintillations. According to these obser- 
vations, for the period of 2001-2005, the slow SW cone was 
~ 56 + 6°, while for the period of 2006-2008, the cone was al- 
most ~ 36 ± 6°. During the intermediate period of 2005-2006, 
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V2 TS crossing on 2007 August 30 
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Fig. 2. Comparison of the V2 plasma and field measurements 
(gray) and MHD model (black solid) for simulation parameters 
a = 42°, /3 = 0°, and ISMF strength 2.2yuG that correspond to 
the best solution obtained in Figure 1 for the Ulysses He param- 
eters. We also show MHD model (black dashed) for simulation 
parameters a = 45°, /3 = 0°, and ISMF strength 3.0 /iG that cor- 
respond to a so lution recently proposed when using IBEX-b ased 
He parameters (McCom as et al. 2012tlBzowski et al. 20121 ). The 
SW speed is ~ 420km/s. From top to bottom, we show the 
plasma speed, density, and magnetic field strength. 
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Fig. 3. (Top). Comparison of the angular positions (ecliptic co- 
ordinates) of computed (solid line) and observed (stars) IBEX 
ribbon estimated for the best fit derived in Figure 1 with SW pa- 
rameters that best correspond to V2 TS-crossing time in 2007 
(SW speed ~ 420km/s). (Bottom). Same but for SW parameters 
that best correspond to the VI TS-crossing time in 2004 (SW 
speed ~ 466km/s). 



the slow S W sector was assumed to shrink linearly between the 
two values (iTokumaru et al. 2010l) . 

In addition to the two values assumed for the slow SW ex- 
tent, we also considered several speeds at 1AU in the range of 
420 to 500 km/s for the slow component in order to cover most 
values thus far measured by the plasma detector on V2 and by 
the Advanced Composition Exoplorer (ACE), an Earth-orbiting 
satellite. For the fast SW co mponent, a factor ~ 1 .9 times the 
slow speed value is assumed (ITokumaru et al. 20101) . For the TS 
crossing time reported so far for VI & V2, the slow solar wind 
plasma density varies in the range of «sw ~ 4.5 - 6.5cirT 3 , 
but was fixed at «sw = 5.2 cm -3 so that the corresponding SW 
dynamic pressure (in the range 1.4 - 2.1 nPa at 1AU) follows 
the square of the SW speed variation. The solar equator is tilted 
7.25° with respect to the ecliptic plane. The ecliptic longitude of 
the Sun's ascending node is 75.77° (J2000). The Parker's spiral 
model, which is assumed for the interplanetary magnetic field, 
has a radial component equal to 35.5 pG at 1 AU. 

In order to compare our MHD simulations to IBEX ribbon 
observations, we use t he locations of enhan ced emission of ENA 
fluxes as reported by dFunsten et al. 20091) . These locations are 
defined by ecliptic latitudes and longitudes at which we com- 
pute the B ■ R = condition of perpendicularity between the 
line of sight of the observation and the magnetic field vector in 



the region outs ide the heliopause where the ISMF has its m axi- 
mum strength (Heerikhu isen et al. 2010tlStrumik et al. 201 It 
is important to stress that modeling of the ENAs flux requires 
kinetic models for neutra ls and is outside the scope of our study 
dHeer ikhuisen et al. 20101) . 

In a first step, we simultaneously target the IBEX observa- 
tion and the V2 TS crossing for which both plasma and field 
measurements exist. The angular grid for ISMF assumes a five 
degree regular step and strengths varying in the set 1.5, 1.8, 2.0, 
2.4, and 3.0 pG. After an approximate solution is derived for 
the strength and orientation of the ISMF vector, we consider a 
finer grid (~2 degrees step) around that solution plus an addi- 
tional value for the strength. Next, we consider the VI TS cross- 
ing in 2004, using those SW speed and density values that are 
consistent with V2 plasma measurements obtained at the time 
of the VI crossing and with ACE measurements conducted at 1 
AU a year or so before its propagation to the VI position. The 
approach described above should provide an ISMF vector con- 
sistent with all observations. 
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3. Results 



Following the method described above, our first step is to 
conduct a sensitivity study by comparing our MHD simula- 
tions with the V2 TS crossing on 2007 August 30 and the 
IBEX ribbon observation. In Figure 1, we compare ecliptic lat- 
itudes/longitudes regions for which the MHD model fits the 
plasma parameters distribution during the V2 TS crossing for 
different strengths of the ISMF. We note that the ensemble 
of solutions is very sensitive to the field strength and orienta- 
tion. However, where the MHD model generally recovers the 
IBEX ribbon, the space regions are not very sensitive to the 
fiel d strength, a finding that is consistent with previou s stud- 



ies ([Strumik et al. 201 lt THeerikhuisen & Pogorel ov 201 lh . It ap- 
pears that while the IBEX ribbon strongly constrains the ISMF 
orientation, it is less effective on the strength, which is better 
constrained by the V2 TS crossing. This means that we need 
both data sets to uniquely derive the ISMF strength and orien- 
tation. This statement is crucial for understanding the difference 
between our results and studies using the only IBEX ribbon data 
analysis. For clarity, a detailled comparison with the recent re- 
sults of Mccomas et al. (2012) is provided at the end of this 
section as a reference. 

Following the proposed step-by step analysis, our parametric 
study shows that our MHD model reproduces the position of the 
IBEX ribbon and the TS crossing conditions for V2 at ~ 83.7 
AU, when the ISMF magnitude is in the range 2.0 - 2.4 fiG and 
is pointing from a direction that corresponds to MHD space pa- 
rameters a ~ 40 - 45° and /3 ~ (-5) - 5°. This solution was 
obtained for a slow solar wind sector of 36° with a density of 
~ 5.2cirT 3 and a speed of ~ 420 km/s at 1 AU, consistent with 
measurements obtained respectively by V2 and by ACE before 
SW propagation to the TS crossing. To conclude this first step, 
we consider a finer grid in the ISMF space parameters and an ad- 
ditional strength of 2.2 fiG as revealed from the previous analy- 
sis. As shown in Figure 1, the best fit is obtained when the IBEX 
solution crosses the V2 TS solution for a ISMF vector with a 
strength of 2.2 + 0.1 fiG and pointing from ecliptic coordinates 
(224,36) + 3°. This solution corresponds to space parameters 
a ~ 42 + 2° and j3 ~ 0.0 + 2° of our MHD simulation. The fit 
to the plasma parameters during the V2 TS crossing is shown 
in Figure 2 (solid line), and the fit corresponding to the IBEX 
ribbon is shown at the top of Figure 3. 

The second step in our approach is to determine whether 
given the SW conditions that best correspond to the TS crossing 
on 2004 December 16, and with the ISMF vector derived in step 
1 , our MHD model provides a good fit to the position of the TS as 
measured by VI in the interplanetary magnetic field. As shown 
in Figure 4 (solid line) a rather good fit is obtained with the same 
ISMF strength and orientation obtained from V2 and IBEX data 
analysis but for SW speed ~ 466 km/s and density 5.2 cur 3 (SW 
ram pressure ~ 1.886 nPa). In addition, for the VI plasma con- 
ditions derived above, a good fit is also obtained for the IBEX 
ribbon (Figure 3, bottom), confirming the weak sensitivity o f 
the outer heliosheath to SW conditions dlzmodenov et al. 2 005). 
However, our sensitivity analysis shows that for the VI events 
thus far reported in mid- 2002 and early 200 3 with the corre- 
sponding noisy IMF data dBurlaga et al. 20031) . the TS position 
could be in the range 85-94 AU. Apparently, predicted posi- 
tions of the TS are at or a few AU ahead of the VI position, 
yet still close to it. Because of large uncertainties in plasma 
and field parameters, our results cannot distinguish between 
the two scenarios proposed so far for t he 2002/2003 event s 
dKrimigis et al. 2003t iJokipii et al. 2004 iKrimigis et al. 2005h . 
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Fig. 4. Comparison of VI field measurements (gray) and MHD 
model (solid line) for simulation parameters a = 42°, /3 = 0°, 
and ISMF strength 2.2 //G that correspond to the best solu- 
tion obtained in Figure 1 but with SW speed ~ 466km/s and 
Ulysses-based He parameters. We also show MHD model (black 
dashed) for simulation parameters a = 45°, [3 = 0°, and ISMF 
stren gth 3.0 uG that corres pond to the IBEX-based He parame- 



ng 

ters dMcComas etal. 20121) . 



Ultimately, detailed modeling using the ISMF vector proposed 
here will be much needed in the future in order to check the flux 
level and energy range of particles measured by V 1 . 

To clarify why it is critical to use most of the available in situ 
measurements and the remote observations to properly constrain 
the ISMF vector, we decided to compare our findings to the re- 
cent report on the ISMF strength and orientation derived from 
the only IBEX data set, including the new ISM He orientation, 
speed, and temperature (hereafter called IBEX-based He param- 
eters) dBzowski et al. 20121: iMcComas et al. 2012b . Our simula- 
tions show that the corresponding MHD model (dashed line) 
does not correctly fit the in situ measurements obtained by V2 
in 2007 (Figure 2) and VI in 2004 (Figure 4). Because the ob- 
tained (dashed) curves fall on opposite sides of the TS posi- 
tion respectively observed by VI and V2, we do not see how 
a stronger ISMF (> 3.0 yuG) would help better fit both Voyager 
measurements. Because the Voyagers' TS crossing and IBEX 
He data were obtained at different epochs, the obtained misfit 
may have several readings. First, our results cast doubt on the 
recently proposed strong ISMF solutions, and the consequent 
conclusions regarding the non-existence of the ISM bow shock 
(IMcComas et al. 2012b . Second, for the newly assumed He pa- 
rameters, other ISMF vector solutions may exist and should be 
tested in the future following the sensitivity approach sketched 
in Figure 1. Finally, if a stronger ISMF is assumed, the mis- 
fit obtained here with the Voyager data could also indicate 
that either the IBEX He data should be critically re-analyzed 
(Mob ius et al. 20121) . or the IBEX He parameters reflect a real 
change in the ISM flow that occurred after the V2 TS crossing 
in 2007. 

To further verify and extend the scope of our results, we 
tried to check their consistency. First, we point out that our 
analysis put the ISMF orientation very close to the observed 
HDP plane as reported by Lallement et al. (2010). Furthermore, 
both the orientation and strength derived here are consistent 
with the values obtained from an independent analysis con- 
strained by the same HDP plane and a full kinetic treatment 
of neutrals dHeerikhuisen & Pogorelov 201 lh . This consistency 
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Galatic Coordinates 





r.3.0.. 

.60... 



Fig. 5. Galactic coordinates of ISMF. (Yellow, hashed square) 
Sky background Lyman-a solutions (37+14°, 49.5+8.5°). (Red, 
small circle) Solutions from Voyager & IBEX data analysis (28, 
52)±3° (this work). (Blue, la rge hashed circle) Starlight polariza- 
tion solutions (38, 23)+35° dFrisch et al. 201 Oh . Big cross shows 
the solution of the Faraday rotation mea sure and dispersion mea- 
sure (5°, 41°) for 7 pulsars dSalvati 20101) . Filled diamond is up- 
wind direction. Angles increase counterclockwise. 



between independent studies justifies our assumption of a con- 
stant H flux. In addition, we tried to compare the ISMF vec- 
tor derived here to existing estimates obtained by other dif- 
ferent techniques. Table 1 summarizes most of the reports on 
the ISMF derived from the VI & V2 TS crossing and IBEX 
ribbon interpretation. Independently, Lyman-ff sky-background 
observation, rotation measure and dispersion measure of ra- 
dio pulsars, along with starlight polarization, provided an es- 
timati on of the ISMF vecto r in the sola r system neighbor - 
hood dRatkiewicz et al. 2008t ISalvati 20 id iFrisch etal. 20101) . 
As shown in Figure 5, the comparison of Galactic coordinates for 
the different estimates shows good convergence between the dif- 
ferent techniques, lending credence to our results and the method 
utilized to obtain them. 



4. Summary and conclusions 

We have shown that the asymmetries observed in the sky 
background Lyman-Q' emission deep in the heliosphere, in the 
TS position as measured by VI and V2, and in the IBEX 
ribbon of ENAs, can all be explained with a rather weak 
|Boo| ~ 2.2 + 0.1 fiG ISMF pointing from Galactic coordinates 
(28,52) + 3°. This direction makes an angle of ~ 42° with the 
Ulysses helium flow in agreement with the value first found by 
dBen-Jaffel et al. 20001). The ISMF vector derived here lies in the 
HDP plane ( Lallement et al. 20101) . consistent with an indepen- 
dent analysis of the SOHO/S WAN hy drogen flow deflection an- 
gle dHeerikhuisen & Pogorelov 201 lh . We have also shown that 
the magnetohydrodynamic modeling of the only V 1 and V2 in 
situ measurements or IBEX data set does not lead to unique val- 
ues of the ISMF strength and orientation as previously claimed 
(e.g., Figure 1). The asymmetry thus far reported for the TS 
shape is not only the result of the ISMF effect but also of the 
SW anisotropy and time variation. 

Our solution for the ISMF vector indicates a weak field that 
is consistent with all reported in situ measurements and remote 
observations. For the Ulysses-based He parameters thus far as- 
sumed for the interstellar flow, a bow shock must exist. By con- 
trast, using a strong magnetic field (> 3.0/iG) making ~ 45° 



from the new IBEX-based He flow leads to a plasma configura- 
tion that is not consistent with the V2 and VI TS crossing. Our 
results cast doubt on the strong ISMF solutions recently pro- 
posed by McComas et al. (2012). Yet for the newly proposed 
interstellar He parameters, further simulations are required be- 
fore one might conclude whether the interstellar bow shock does 
disappear under those assumptions. 

The ISMF derived here, at the actual Galactic position of the 
Sun, is in the range ~ 70 - 80° from the different estimations 
of the large-scale (kpc) Galactic magnetic field (e.g., Figure 5). 
This finding is consistent with previous studies of rotation mea- 
sure and dispe rsion measure of radio pulsars and from starlight 
polarization dSalvati 20101 : IFrisch et al. 20101) . 
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